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SOME LIMITATIONS ON THE USE OF DAMPING IN SHORT PRESSURE PROBES 
by Ted W. Nyland and Robert C. Anderson 
Lewis Research Center 

SUMMARY 

An analytical study is described in which variations in the useful frequency response 
of a short pressure probe are calculated as a function of viscous losses (damping) and 
operating conditions. The purpose of the damping is to increase the useful frequency 
response of the probe where useful response is defined as that frequency range over 
which the amplitude ratio is 1.0±0.05. The study is based on the use of the Bergh and 
Tijdeman recursion formula, which was modified to allow for the inclusion of sections 
with sets of parallel tubes. A computer program used to calculate the response of 
probes with the modified recursion formula is given in the appendix. Results indicate 
that damping can be used to increase the useful response of a probe by a factor of three 
or more. The gain in useful response, however, is highly dependent on the operating 
pressure and temperature level. For a temperature of 300 K and average pressure of 

10 newtons per square centimeter, a factor of three in useful response is achieved only 
over a 15-percent variation in the average pressure. At 300 K, the size of capillary 
tubes required to get useful damping is strongly dependent on pressure level. Tube 
sizes used in this study ranged from an 0. 013 -centimeter inside diameter for use at 

11 newtons per square centimeter to 0.0088-centimeter inside diameter at 26 newtons 
per square centimeter . 


INTRODUCTION 

A number of experiments being conducted at the Lewis Research Center have re- 
quired the measurement of transient pressures. Typical of these experiments is one in 
which the performance and stall limits of a turbofan engine (ref. 1) were determined. 

In this experiment, transient pressures were measured with probes consisting of short 
lengths of small-diameter tubing connected to a housing that held a miniature pressure 
transducer, resulting in a lightly damped measuring system. All of the probes in this 
experiment were of a minimum length as determined by the physical constraints imposed 



by the test environment and were designed for maximizing the natural frequency. For 
this experiment, it had been estimated that a useful response of 500 hertz was required, 
where useful response is arbitrarily chosen to be that frequency at which the amplitude 
ratio is 1.00 within ±5 percent. In other test programs being conducted at Lewis, use- 
ful frequency responses up to 1500 hertz are desired. To extend the useful response of 
many of these probes, a potentially attractive method involves increasing the viscous 
damping in the probe. 

Generally speaking, most transient pressure probes can be considered to be lightly 
damped, second-order systems. Typical probes have undamped natural frequencies in 
the 2000- to 3000-hertz range with useful frequencies of 20 percent of the natural fre- 
quency (see ref. 2 for second-order response equation). With optimum damping, the 
useful frequency of a second-order system can be increased to better than 80 percent of 
the natural frequency. This more than quadruples the useful response range over that 
of an undamped probe and thus offers an attractive solution to the problem. It should be 
pointed out that the increase in useful response is not achieved without some sacrifice. 
The sacrifice in this instance is that a phase shift is introduced into the measurement 
which may or may not be acceptable in given measurement situations. 

Experimental data have shown that the simple second-order system equations do not 
accurately predict the amplitude response over the entire frequency range of probes of 
complex geometry and/or those with appreciable damping. However, Iberall (ref. 3) 
developed an equation, not based on a second-order system approach, from which the 
amplitude response of a single tube and volume pressure transmission system can accur- 
ately be calculated. Iberall' s equation is derived from the fundamental flow equations 
(that is, the Navier-Stokes equation, the equation of continuity, the equation of state, and 
the energy equation) and assumes that the sinusoidal disturbances are very small com- 
pared with the static conditions in the system. Goldschmied (ref. 4) has shown that 
available data correlate well with Iberall's equation. Following Iberall's approach, 

Bergh and Tijdeman (ref. 5) have developed a recursion formula from which the response 
of a series connected set of tube and volume systems can be computed. The Bergh and 
Tijdeman formula reduces to the Iberall equation for the case of a single tube and volume 
system. References 5 and 6 indicate that good agreement is achieved between experi- 
mental data and the recursion formula for many types of tube systems. As the recursion 
formula appears to adequately predict the response of many types of pressure measur- 
ing systems, it, along with a modified form discussed later, was used for all computa- 
tions described herein. 

An analytical study was made at Lewis for the purpose of determining what gains 
could be made in the useful response by increasing the damping in short pressure probes 
and examining some of the problems that might arise in practice. It was also of interest 
in the study to determine what changes would occur in the useful frequency response when 
operating at off -design temperature or pressure conditions in a probe. This report is a 
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summary of that study. A discussion of some of the practical problems associated with 
the use of various damping materials is also presented. 




ANALYTICAL MODELS AND RESPONSE EQUATIONS 

t, A diagram of the basic probe used in the study is shown in figure 1. This probe 

consists of a 2. 54-centimeter long section of tubing connected to a small cavity contain- 
i ing a miniature pressure transducer. The inside tube diameter is 0.325 centimeter. 



Figure 1. - Diagram of basic probe. (Dimensions are in centimeters. I 


The volume ratio of the probe, defined as the ratio of cavity volume to tube volume, is 
0.027 and is considered representative of a minimum volume ratio that could be achieved 
in an actual probe design. The amplitude and phase angle response of this probe were 
computed from the recursion formulas of reference 5. The subroutines used in these 
computations are listed in reference 6. 

The useful response of the basic probe is shown in figure 2. Plotted are curves of 
constant useful frequency as a function of pressure and temperature. As an example, at 
300 K and 10 newtons per square centimeter, the useful response is seen to be 
625 hertz. Any values of useful response discussed in this report can be referenced 
to the values in this figure since these represent the response of an undamped probe. 

Two approaches could conceivably be used to increase the damping in the basic 
probe. In the first approach, one might consider simply working with smaller diameter 
probes and determining their effect on the useful response. The second approach would 
involve filling the probe tube with an insert of some material that would increase the 
viscous losses while maintaining a large cross sectional area for flow. The insert might 
take the form of a porous type of material such as compressed steel wool or it might 
take the form of a set of small-diameter parallel tubes. In both forms, the flow area 
would be made up of a large number of small-diameter passages that would result in 
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Figure 2. - Lines of constant useful frequency response 
as a function of temperature and pressure of an un- 
damped probe. 
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Figure 3. - Basic probe with parallel-tube insert section. (Dimensions are in 
centimeters. ) 


high viscous energy losses. For this report, three parallel-tube inserts were examined; 
their diameters, lengths, and numbers are given in figure 3. These inserts replace the 
2.54-centimeter-long by 0.325-centimeter-diameter tube in the basic probe. 

For computational purposes, the flow characteristics of the parallel- tube insert can 
be introduced into the recursion equations of reference 5 as follows. The general sys- 
tem, consisting of N tube and volume sections, used in the derivation of the recursion 
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Figure 4. - Tube-volume system of Bergh and Tijdeman derivation. 
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formula is shown in figure 4. In addition to the variables used in the derivation, we 

define a variable C., which is equal to the number of parallel tubes associated with the 
th J 

j section of the system. Likewise, there is a variable for the j+1 section. 

These variables are included in the derivation by simply multiplying equations (42) 

and (45) of reference 5 by C, and C. 1 , respectively. Equation (42) represents the 

^ th 

total mass flow leaving the tube set in the j section; equation (45) represents the mass 
flow entering the tube set in the j+1 section. These variables are then carried through 
the remainder of the derivation, producing a modified form of the recursion formula: 


J_ = 


P. 


j-1 


V 

v i 
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(The symbols are defined in the appendix. ) A set of FORTRAN IV subroutines based 
on the preceding modified version of the recursion formula is given in the appendix. 

For all computations, except those discussed in the section Analysis of a Typical 
Damped Probe Application, N was equal to 1. In computing the frequency response of a 
given probe, the density of air at various temperatures and pressures was calculated 
from the perfect gas law equation. The air viscosity was calculated from an equation 
given in table 1-b of reference 7 (p. 10). A summary of these equations and other 
property values used in the recursion formulas are given in table I. 


TABLE I. - PROPERTY VALUES AND EQUATIONS 


USED IN COMPUTATIONS 


Temperature, K 

T 

2 

Pressure, N cm 

P 

Specific heat ratio 

1.4 

Prandtl number 

0.7 

3 

Density, g/cm 

3.485X10' 2 P 'T 

Viscosity a , g/(sec)(cm) 

1.458X10' 5 T 1 ' 5 /(T + 110.4) 


a From ref. 7. 


EFFECTS OF TUBE DIAMETER ON THE USEFUL FREQUENCY OF PROBES 

On first thought, it would seem reasonable that the useful frequency response range 
of a transient pressure probe could be increased by using a smaller diameter tube in the 
probe. The small tube would result in an increase in viscous energy losses and there- 
fore an increase in the effective damping of the probe. However, to gain in useful re- 
sponse , one must maintain a high natural frequency . The natural frequency of a probe 
is dependent on volume ratio and probe length. 
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As indicated in reference 6, to maintain a high natural frequency for a given length 
probe requires that the cavity volume be reduced proportionately with reductions in the 
square of the tube diameter. If the volume ratio is held constant, an increase in the 
useful response range can be achieved using smaller diameter tubes. This is illus- 
trated in figure 5, which shows curves of the calculated useful frequency as a function of 
tube diameter for constant volume ratio and for constant transducer cavity volume . As 
is shown by the curve for constant volume ratio, the useful frequency of the basic probe 
could be optimized by using tube diameters of about 0.012 centimeter while maintaining 
a volume ratio of 0.027. However, it should be pointed out that this curve is represen- 
tative of an ideal probe design, as one cannot build such a probe when using available 
miniature pressure transducers and allowing reasonable machining and construction 
tolerances. 

A more realistic design condition is represented by the second curve in figure 5. 

This curve shows the computed useful frequency for the basic probe assuming that the 
volume is held constant while the tube diameter is varied. It is seen that the optimum 
response is achieved when using as large a diameter tube as possible and that no increase 
in the useful frequency can be obtained when using small diameter tubes. 



Tube diameter, cm 


Figure 5. - Useful frequency as function of tube diameter. 
Temperature, 295 K; pressure, 9.92 newtons per square 
centimeter. 
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FREQUENCY RESPONSE CALCULATIONS FOR PROBES 
WITH PARALLEL -TUBE INSERTS 


Typical response curves calculated with the modified recursion formulas for a probe 
with aparallel-tube insert are shown in figure 6. Plotted are the amplitude and phase 
angle response curves for the basic probe and for the basic probe with the 1-204 and 
1-639 inserts. The total cross sectional area for flow of the basic probe with and without 
the inserts are equal. The experimental data points are for the basic probe taken from 
reference 6 and show good agreement with the computed curve. As would be expected, 
increasing the number of tubes while decreasing their diameter results in a lowering of 
the peak amplitude ratio. In general, it is seen that, by use of different inserts, the 
amplitude response can be made to vary in a manner similar to that which occurs when 
the damping coefficient in a second-order system is changed. Both a reduction in the 
amplitude ratio and the introduction of a phase shift occurs. For the basic probe and the 
1-204 insert, the intersection of the amplitude curves with the ±5 percent (±0.42 dB) lines 
show useful frequency responses to about 600 hertz. By using the 1-639 insert, the use- 
ful response has been extended to about 1300 hertz and results in more than doubling the 



Figure 6. - Frequency response curves for basic probe and probes 
with 1-204 and 1-639 inserts at temperature of 300 K and pres- 
sure of 10 newtons per square centimeter. 
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response range of the basic probe. Thus, it appears that the useful response of a probe 
can be increased by use of these inserts. 


TEMPERATURE AND PRESSURE EFFECTS ON THE 
RESPONSE OF A DAMPED PROBE 

The amount of damping achieved with a given insert is a function of the density and 
viscosity of the fluid in the probe. Since these properties depend on temperature and 
pressure, the useful response of the probe also is temperature and pressure dependent. 
Figure 7 shows some typical response curves for a probe with an 1-639 insert. In fig- 
ure 7(a), response curves for the probe operating at 250, 300, and 350 K are shown for an 
assumed constant pressure of 10 newtons per square centimeter. The useful response 
is seen to change from 1000 to 1275 to 750 hertz at these temperatures. In figure 7(b), 
response curves for pressures of 5, 10, and 15 newtons per square centimeters are shown 
for an assumed constant temperature of 300 K. For these pressures, the useful re- 
sponse varies from 300 to 1275 to 750 hertz. The useful response for the basic probe is 
around 650 hertz. Thus, there exists a range of limited temperatures and pressures 
over which an insert such as an 1-639 could be used for extending the useful frequency 
range of a transient pressure measuring probe. 

A more complete picture of the range of temperatures and pressures over which a 
probe with an 1-639 insert could find use is shown in figure 8. Plotted are lines of con- 
stant useful frequency as a function of temperature and pressure. As is shown in the 
figure, there is a region in the temperature and pressure domain in which the probe’s 
useful response is greatly increased. As the operating temperature and pressures in- 
crease, the range in temperature and pressure over which this probe could be used also 
increases. However, at temperatures such as 300 K, the probe is limited to operating 
at 10±1. 5 newtons per square centimeter while having a useful response of 1800 hertz or 
more. For general probe applications, this small allowable variation in operating pres- 
sure imposes a severe limitation on the usefulness of a damped probe. 

Figure 9 is a similar plot of lines of useful frequency as a function of pressure 
and temperature for a probe with an 1-1367 insert. The overall shapes of the lines 
for this probe are the same as those shown in figure 8 but are more broadly spaced 
and displaced to the right along the pressure axis. At 300 K, this probe could find use at 
25±3.0 newtons per square centimeter and still maintain a useful response to at least 
1800 hertz. Again, this is a fairly small allowable operating pressure range and se- 
verely restricts the use of damped probes of this type when operating at low tempera- 
tures. Even with a relaxation of the useful response to 1200 hertz (approximately doub- 
ling the range of the basic probe), the pressure range for operation would be increased 


9 



/-+5 percent 



J I I I 1 I I ll I I I I ! I I ll 





Frequency, Flz 


(bl Showing effects of pressure. 

Figure 7. - Frequency response curve for a probe with an 1-639 
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Figure 8. - Lines of constant useful frequency as func- 
tion of pressure and temperature for a probe with an 
1-369 insert. 


Figure 9. - Lines of constant useful frequency as function 
of pressure and temperature for a probe with an 1-1367 
insert. 


only to ±2 newtons per square centimeter for the 1-639 insert and to ±6 newtons per 
square centimeter for the 1-1367 insert. In general then, for measurements around 
300 K, it does not appear practical to use damping to increase the useful response of 
probes unless the average pressure is known and tightly held to within 10 or 20 percent 
of the probe design value. At higher temperatures, this restriction in operating condi- 
tions becomes less severe and might allow one to use damping for increased frequency 
response. We feel that any type of materials used to increase the damping in a probe 
will have similar operating restrictions and that, in general, these results presented 
for the parallel-tube inserts are applicable to a wide variety of damping materials. 


ANALYSIS OF A TYPICAL DAMPED PROBE APPLICATION 

In turbojet engine experiments, transient pressure probes are sometimes positioned 
between compressor stages in regions of elevated temperature. These temperatures 
are generally higher than the allowable operating temperature for the miniature pressure 
transducer, and water cooling is normally used to maintain a temperature that is com- 
patible with the transducer's operating requirements. This temperature is normally 
around 300 K. The temperature at the probe inlet may range from 300 K to more than 
700 K, depending on the probe's position in the compressor. The useful response then 
becomes a function of the temperature distribution between the inlet and cavity volume 
of the probe. One, therefore, must take this temperature distribution into consideration 
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when determining the useful response range of a probe. 

A set of calculations were made to determine what useful response could typically 
be achieved under test condition which resulted in a thermal gradient between the inlet 
and cavity volume of a damped probe. A linear temperature distribution was assumed to 
exist over the probe length with the cavity volume at 300 K and the probe inlet at a tem- 
perature corresponding to the pressure levels as shown in figure 10. The curve in fig- 
ure 10 corresponds approximately to the temperature and pressure distribution that 
would occur in a compressor with an inlet temperature of 300 K, an inlet pressure of 
5 newtons per square centimeter, and an efficiency rj of 88 percent. The results of the 
response calculations are shown in figure 11 where it was assumed, for computational 
purposes, that the probe consisted of 10 equal length sections (N = 10). Plotted are the 
ratios of useful response for probes with the 1-639 and 1-1367 inserts to the useful re- 
sponse of a basic probe. From the curves in figure 11, the useful response of a probe 
operating at average pressures of 20±3 newtons per square centimeter could be doubled 
by use of an 1-639 insert. Similarly, the use of an 1-1367 insert would double the useful 
response range of a probe designed for average pressures of 60±15 newtons per square 
centimeter . These operating levels are somewhat higher than those that would exist if 
the probes were held at constant temperature (see figs. 8 and 9). This points out the 
necessity of knowing the temperature distribution in a probe when one is considering the 
application of damping to short probes. As was the case for the constant -temperature 
probe at low pressures, the allowable operating range is very small and may limit the 
general application of the inserts. At high pressures, one could consider the use of in- 



Figure 10. - Probe inlet temperature as function of 
average pressure level. Compressor inlet tempera- 
ture T c , 300 K; compressor inlet pressure P c , 5 
newtons per square centimeter; compressor efficiency 
H, 0.88; specific heat ratio y, 1.4. 



Figure 11. - Example of change in useful response as 
function of pressure level for two damped probes as- 
suming a linear temperature gradient between probe 
inlet and cavity volume. Probe inlet and cavity vol- 
ume temperatures are given in figure 10. The re- 
sponse ratio is defined as the ratio of useful fre- 
quency for a damped probe to that of the basic probe. 
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serts to increase the useful response by at least a factor of two over that obtained with 
a probe with no insert. It is interesting to note that by use of an 1-639 insert, the useful 
response would be equal to or higher than that achieved with an undamped probe at all 
pressures above 12 newtons per square centimeter over the range indicated. The same 
is true for an 1-1367 insert above 35 newtons per square centimeter. 


INSERT DESIGN CONSIDERATIONS 

Preliminary screening of materials that may be suitable for use as inserts can be 
made by considering two basic design requirements. The first requirement is that the 
probe and insert combination should have a small volume ratio. For a damped probe, 
the volume ratio is defined as the ratio of cavity volume to the total volume of the flow 
passages in the insert. The reason for this requirement is illustrated in figure 12. 



Figure 12. - Useful frequency as function of number of 
OL 01286-centimeter-inside-diameter tubes at 330 K and 
10 newtons per square centimeters. 


Shown is a plot of the useful response of a damped probe as a function of the number of 
0. 0128-centimeter -diameter tubes in the insert. Also indicated along the curve are cor- 
responding values of volume ratio. (For these calculations, the volume ratio is in- 
versely proportional to the number of tubes.) Initially the curve rises sharply and be- 
comes asymptotic to some value of frequency at which large increases in the number of 
tubes (decreases in the volume ratio) result in only small increases in the useful re- 
sponse . Depending on how much reduction in the optimum useful frequency response is 
acceptable, one can place lower limits on the number of tubes that would be required in 
the case of a parallel-tube insert, or upper limits on the volume ratio of the probes with 
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porous material inserts. For most damped probe designs, it is recommended that one 
consider only materials that, when used in a probe, would result in volume ratios in 
the 0. 05 range. This 0. 05 volume ratio would result in a decrease in response of not 
more than about 15 percent compared with a system with zero volume ratio. 

The second design requirement is related to the fact that the individual passage- 
ways making up the total flow area be small enough to obtain adequate viscous losses 
in the system. For a parallel-tube insert, the required tube diameters fall into the 
0. 01 -centimeter -diameter range, their actual size being dependent on the operating tem- 
perature and pressure. As shown in figures 8 and 9 a design temperature of 300 K and a 
pressure of 11 newtons per square centimeter would require tubes of 0. 01 2 8 -centimeter 
inside diameter. Changing the design pressure to 26 newtons per square centimeter 
would require tubes of 0. 00879 centimeter diameters. These tubes are small enough so 
that manufacturing tolerances become a problem. 

Other than an actual parallel -tube set, inserts could be made from sintered metal - 
fiber materials. The actual flow passages in these materials are irregular, consisting 
of spaces between the metal fibers. The hole sizes in these materials would generally 
have to be smaller than those associated with a parallel -tube insert in order to achieve 
similar damping effects. So, in general, one can consider using materials that have 
large cross sectional areas for flow and with hole sizes in the range of 0.01 centimeter 
for a tube system and less than that value for a fiber metal type material. 

From an analytical standpoint, it is very difficult and beyond the scope of this report 
to introduce the damping characteristics of porous materials into the recursion formula. 
Any flow parameters specified by the manufacturer of porous materials are normally 
measured under steady-state flow conditions. These parameters cannot readily be con- 
verted to useful parameters for the case of an oscillatory flow. One, therefore, would 
probably have to measure the frequency response of the probe with a particular insert. 
Then by trial and error, the modified recursion equation could be used to determine an 
equivalent parallel-tube insert whose frequency response curve matched the experimen- 
tal data points. Then, trends in the response characteristics for other temperature and 
pressure conditions could be obtained through calculations based on the parallel-tube 
equivalent. 


CONCLUDING REMARKS 

Ideally, the useful frequency response range of a transient pressure measuring 
probe can be increased by introducing damping into the probe body. For a given size 
probe, it has been shown that just decreasing the probe tube diameter will not result in 
an increase in the useful response. At small tube diameters, the volume ratios becomes 
the controlling parameter and, therefore, the cavity volume must be decreased propor- 
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tionally with the change in tube diameter to achieve a gain in the useful response. Inserts 
consisting of sets of parallel tubes can be used to increase the viscous losses in probes. 
For measurements at standard conditions, diameters for these parallel tubes would be 
in the 0. 01 -centimeter range. The useful response of probes with inserts becomes 
strongly dependent on the operating temperature and pressure of the gas in the probe. 

At 300 K and 10 newtons per square centimeter, one can raise the useful response by a 
factor of three if the operating pressure is held within 15 percent of the design pressure. 
At higher temperature and pressure levels, this operating restriction becomes some- 
what less severe. For probes with inserts that are not of the parallel -tube type, the 
frequency response will normally have to be determined experimentally. Trends in the 
useful response curve can then be predicted by use of the recursion formula. When 
damping is introduced into probes, one must remember that a phase shift is introduced 
into the measurement that may or may not be acceptable in all measurement situations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 28, 1971, 

126-61. 
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APPENDIX - COMPUTER SUBROUTINES FOR CALCULATING THE RESPONSE 
AT A FREQUENCY OF A MULTIPLE TUBE-VOLUME SYSTEM 
INCLUDING PARALLEL TUBE SECTIONS 

A set of subroutines for computing the amplitude ratio and phase angle at a fre- 
quency for a multiple tube-volume system with parallel tube sections is described 
herein. The subroutines are based on the modified Bergh and Tijdeman equation as 
discussed in a previous section of this report. The symbols used in the modified equa- 
tion and computer program are defined as follows: 


Engineering 

Computer 

Definition 

symbol 

symbol 


— 

A 

amplitude ratio 

C 

CONGL 

number of tubes in given section 

D 

D 

tube diameter 

f 

FREQ 

frequency 

i 



J n 


Bessel function of first kind of order n 

L 

XL 

tube length 

N 

N 

number of tube -volume sections 

P 

— 

amplitude of pressure disturbance 

Pr 

P 

Prandtl number (usually equal to 0.7) 

— 

PA 

phase angle 

P 

s 

AMPRES 

mean pressure 

V t 


volume of a single tube in a given section 

V v 

V 

cavity volume of given section 

y 

GAMMA 

specific heat ratio 

n 

VIS 

fluid dynamic viscosity 

Ps 

RHO 

mean mass density 

(g) 

DVDP 

change in volume V due to change in pres- 
sure 
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u)=2ni — angular frequency 

j I tube -volume section under consideration 

The FORTRAN IV subroutines are referenced by the following call statement: 

CALLAMPPHS (FREQ, A, PA, I, N) 

where FREQ is the oscillatory frequency for which the amplitude ratio and phase angle 
are to be computed, N is the total number of series connected tube -volume sections, 
and I is the section for which the amplitude ratio and phase angle are being computed. 
FREQ, N, and I must be defined before calling the subroutine. The computed ampli- 
tude ratio P(I)/P(0) will be placed in A, and the computed phase angle in PA on com- 
pletion of a computation. 

The program accepts as many as 10 series -connected tube -volume sections in a 
system. The parameters for each section must be put into the subroutines through three 
common blocks referenced as follows: 

COMMON/GEOMTY /D(10) , XL(10), V (10), DVDP (10) 
COMMON/THERMO /GAMMA, P (10), RHO (10), VIS (10), AMPRES 
COMMON/CON/CONGL (10) 

The geometry parameters, diameter (D), length (XL), volume (V), and volumetric dis- 
placement (DVDP), for each section are entered through the GEOMTY common block. 

Zero values (0.0) may be entered for any or all values of V and DVDP. The thermo- 
dynamic fluid properties, ratio of specific heat (GAMMA), density (RHO), viscosity (VIS), 
and ambient pressure (AMPRES), are entered through the THERMO common block. The 
ratio of specific heats and ambient pressure are assumed to be equal in all sections. The 
number of parallel tubes (CONGL) in each section is entered through the CON common 
block. For a single tube section, a 1.0 must be specified for the CONGL value asso- 
ciated with that section. 

The program used for generating the value of the Bessel function is limited in use 
to Bessel function of integer order with arguments of the form C (1 - i) where C is a 
real constant. Note: In computing the value for a Bessel function, very small numbers 
are generated that may result in machine underflows . It is suggested that use be made 
of a major underflow option if available on the machine being used. These underflows 
have been found to occur in the evaluation of statements 1000 and 2000 in the SERIES 
subroutine . 
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A listing of the nine subroutines required for computing the amplitude ratio and phase 
angle at one frequency for a series connected system of tube- volume sections including 
sections with parallel tube sets using the modified Bergh and Tijdeman equation is given 
next. 


JIB FTC SUE A m p LIST 

SUER (JUT I.M^ AMPPhS! FULCRA, PA, I, J) 

C 

C 

C THIS ROUTINE CALCULATES THE AMPLITUDE RATIO (A) P ( I > / P C 0 ) 

C A'Jf,' PHASF A ‘ i 0 L F (PA) FCR A SERIES CONDUCT I CN OF J TUBE-VOLUML 
C SYSTEMS GIVEN F R E OU F A C Y ‘FREQ*. 

C THE SYSTEM GECMFTRY IS ENTERED THROUGH COMMON BLOCK /GFOMTY/. 

C THF THERMODYNAMIC CONDITIONS ARE ENTERED THROUGH COMMON BLOCK / THFRMO/ . 
C 

c 

DIMENSION PRATIOdC) 

COMMON/GEOMT Y/D ( 1C ) ,XL ( 10) ,V ( 10) ,DVI)P( 10) 

COMMCM/ ThFRMO/GAMM A ,p ( !0) ,RH0( 10 ) , VIS ( 10) tAMPRf.S 

COMMOD/XFREO/FPFCl 

COMMON /CGN/CCNGL ( 10 ) 

COMPLEX PRAT 10, RPCLRS 
FREQ1 = FRFO 
RECURS = 1.0 
DO 10 M = 1 , 1C 
1C PRATIC(M) = C.O 

N = J 

CALCULATE SUCCESSIVE PRESSURE RATIOS 
DO ?C K=l,.v 

CALL RATIO! PRAT IC»J+1-K,J) 

CALCULATE R FLU I R FI) RATI I , HY MULTIPLYING APPROPRIATE PRESSURE RATIGS 
D030 K = 1,1 

RECURS = RECURS* FRA TI C ( K ) 

A = CAPS (RECURS) 

P A= ATA\2 ( AI 'MG ( RECURS ) ,RFAL ( RECURS) ) 

CONVERT RADIANS TO DEGR'ES 

PA= PA* c > 7. 2957705 
IF(PA.GT.C.C) PA = -360. +PA 
R F T U R N 
END 
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TIBFTC XRATIU LIST 

SUt-iT OUT INF PATIf!(PR/;TIO, J,K) 


THIS PROGRAP CALCULATFS THE PRESSURE RATIO P(J)/P(J-1). 


D I w E (\ S I f. i\ PRAT I G ( 10 ) , VP AT I U ( 10 ) , CF’RRCTI 1U ) 

COMUOM/OtGMTY/DI 10) ,XL ( 10) ,V(10) ,0VDP(10) 

C 0 R M 0 M / T h E R p 0 / G A R A ,P { 1C ) , RHO (10) » V I S ( 10) , A P P R E S 
C ON MON/ X FR EC / FR E 0 
COPMCN/CQ M/CCMGL ( 10 ) 

CO’PL Fa PI i I , JR AT , POLY , PCLRAT.CCOSH ,CS IMH, PRAT IC, TERM1 ,TERP2 .TERM 3 , 
CTPHI 1 , TPH I 2 t JRAT If 1 
PI = T. 1415 927 

VRATIf(J) = V( J ) /(P I*D( J)**2*XL( J ) /4. ) 

VRATIO(J) = V R A T I 0 ( J ) / C F N G L ( J ) 

RFRFO = 2.*PI*FREC 
IF (V(J ) , !n E . C . 0 ) G L T 0 1 
CUPRCT ( J ) = 1.0 

G0T03 

1 CI1RRCT ( .1 ) = 1.0 + (LVI)P(J)*AMPRES/V(J) ) 

3 CCMTINUF 

I F ( J . rg.K ) GO TU 5 

COMP AT = r,G\'Gl ( J + l ) /CONOL ( J ) 

POLRAT = CSQR T( PLLY ( .1 ) / POL Y ( J + l ) ) 

RHORAT = SCR I ( R H L ( J + 1 )/RUG(J) ) 

VOL RAT = 0(J+1)*»2/(D(J)**2) 

TPHI? = PHI (J + l ) 

JRAT = CSQPTI 

JJR AT I P ( CSCRT (CMPLX ( C. 0 , - ( REREQ*f)( J + 1 )**2*RHC ( J + l ) / 

J( 4.G*VIS( 1 + 1) ) ) ) ) ) / 

J J R A T I L l C S 0 R T ( CM P L X ( C . 0 , - ( K F R E Q * 0 ( J ) * * 2 * R H C ( J )/ 

J( 4 . 0 * V 1 S ( J ))))))) 

5 TPHI 1 = PhHJ) 

T E R V 1 = CCUSHtTPFI 1*XL ( J) ) 

TERM 2 = V p AT I 0 ( J ) «C PR R C T (J)*XL(J)*TPHI1*PULY(J)*CSINH(TPHI1*XL(J) ) 
1 F { J.Pii.K ) GC TO 2 

TLRM3 - ( CS I N H ( TPh I I*XL ( J ) ) /CS I NH ( TPH I2*XL ( J+l ) ) ) * 

T VOLRAT * R H C R A 1 *JRAT * P 0 L R A T * 

t i consul r PF I? *XL ( J+ 1 ) l-PRAT IOIJ + 1 ) ) 

TERM 3 = T t R M 3 c C 0 N R A T 
GUT 0 6 

2 TCRY3 = C^PLXIO. 0,0.0 ) 

6 CnuTI’.'UF 

PRATIF.(J) = 1.0/ ( TERM1 + T C RM2 + TERK3 ) 

RETURN 

ENL. 
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$ I P F T C r FF 1 F LIST 

SU p R T-U T I IE SEP I FS 
C 

c 

C THIS RCUTINF C/'LCUL A TFS TH r VALUE OF THE SFRIFS N E E 0 E D FOP 
C HANK ELS ASYMPTOTIC EXPANSION 
C 

C IN CALCULATING p ( J ) ANO 0 ( J ) NUMBERS WITH EXPONENTS LESS THAN -38 
C PAY CCCIR. THIS PROGRAM USES A D UII.T-IM P.OliTP'E V H I CM SFTS P C J ) 

C ANC/CR C(J) EQUAL TO Z r RP WUFMFVFR THESE SMAlt NIP'RFRS HCCUR. 

C 

C 

C 

C 

C n M F L F X Z, MU,PSUM, OSUM, FM, P ( 50) , Q( SO ) 

C3NMGN/S 1/X , Y, PSU M , QSIIM 
CANNON /FMN AM/FM 

psl'y = ( o. , ) 

CS UN = ( 0 . , 0 . ) 

YU =4.*FM**? 

7 = C y P L X ( X f Y ) 

CM ICO J=1 , SO 
N = J-l 
FM 

I r ( N .NE. 0) GO TO 30 

20 pi J) =( i.,n. > 

C(J) = ( MU- 1 . ) / ( 3. *7 ) 

GO TO SO 

30 F" 2 = 2.*EN 

F M 4 = 2 . * F N 2 
FN21 = FN2 + 1. 

F 14 1 = 2 . * F N 2 1 

IT (CAPS ( P( J-l )). PO.O. ) GO TO SO 

10C0 P(J) = -(( Mil- ( FN4-3 .)**?)/< 44.* ( FN2-1 . >*FN? ))*(( MU- ( FN4-1 .)**?) /Z 
t**2)*P(J-l) 

If (J .GT. 2 .ANO. C A B S ( P ( J ) ) .GT. CARS(P(J-1>) ) P ( .1 ) = (0.,0.) 
GO TO c 5 

so p ( j ) = (c.,o.) 

55 IT (CAPC (0( J-I ) ) .FP.O. ) on TO OP 

20C r C(J) = — ( ( Ml! — ( FYJ41 — 3. )**2)/(64.*F'|2*FN21) )*( ( M U- ( F N 4 1 — 1 . )**?) / 7 
$ * * 2 ) *C (J-l ) 

IF (J .GT. ? .ANn. CA«S(0(J)) .GT. CAPS(OU-l)) ) 0 ( J ) = (0.,0.) 
GO T 0 0 0 

8 0 C ( J ) = ( 0. , C . ) 
qQ PCL'N = P(J)4PS')M 

osi'N = cmj ) +QSUM 

IF (f APS ( P( J ! ) .FA.'-, .amp. CAPS(Q(J)> .EO. 0.) GO TM no 
100 CONTINUE 
11C CONTINUE 
R n Tl RN 
E‘T 
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tiBf-rc bessnn list 

SUBROUTINE Z ECE E S < J ,N , Z > 


THIS ROUTINE CALCULATES THF- VALUE OF THE BESSEL FUNCTION OF THE 
F I p ST KIND C1F INTEGER 0 P D F R . THE ARGUMENT IS OF THE FORM 
C(1 - SORT ( - 1 ) ) WHERE C. IS A REAL CONSTANT LESS THAN 30 . 

J... BESSEL FUNCTION VALUE 
N.... BESSEL FUNCTION ORDER 
Z... BESSEL FUNCTION ARGUMENT 


COMPLEX J,Z 

D0U3LF PRECISION CF i\, T , R I 2 ) ,C , DK , U , B 

0 = REAL(Z) 

C = D * L .41421 3562 3730950 0 
CR = D 

DFN = F L 0 A T ( N ) 

IF(N.Nc.i) GO TO 1 
T = 5.0-1 
GO TO 2 

1 T = l.DO 

2 R ( 1 ) = T 

R ( 2) = 0.00 

C = ( C * * i . ) / 4. CO 

UK = C.CO 

1 = 1 
M = 0 

5 OK = UK + l.DO 

T = T * C / ( ( UK ) * ( I'.K + DFN ) ) 

1 = 1 + 1 

IF (T.LT.UAB, S(R( I ) * i . 0 - P ) ) GO TO 7 
M = M + 1 

B = 1 . DC 

IF(M.GT.l) B = - 1 .DC 
I F ( M . F 0 . 3 ) M = -1 
R ( I ) = R ( I ) + R * T 

IF! I .FQ.,i) I = 0 
GO TO 'j 

7 RR= R ( 1 ) 

R I = R( 2) 

IF (N.Nl .2) GO TO E 
R I = ( - 1 . ) * R(l) * C 
R R = R ( ^ ) * C 

8 J = C V ‘PLX(RR,RI ) 

RETURN 

END 


STBFTC XJRATu LIST 

COMP L F X (-UNCTION JRATIU(CARG) 

c 

c 

C THIS ROUT IN L CALCULATES THE RATIO OF THE BESSEL FUNCTION OF THE 
C FIRST KINC, ORDER 2 TO 1HC PESS^L FUNCTION OF THE FIRST KIND, 

C OR l.i FR 0 WIlH COMPLFX ARGUMENT. 

C 

C 

C OMP L r X EM, Q SUM, P SUM, XP SUM ( 2) , XQSUM ( 2 ) • XJ( 2) , RATIO, XX, Y 
IV, Z 

COMPLEX J?,J0, JRATIC.CARG 
COMMON /Sl/X, Y, PSUMjOSUM 
COMMON /FMNAM/FM 

DATA PI/3. 141^92?/»M/-l/,Z/(0. 0,1.0)/ 

IFIRFAL (CAR.,) .GT. 30.0) GO TO 40 
C 

C IF T HP REAL PART OF TEE ARGUMENT IS GREATER THAN 30 USE IS MADE OF 
C EA'-.KELS ASYMPTOTIC EXPANSIONS (SEE NBS HANDBOOK OF MATHEMATICAL 
C FUNCTIONS AMS M!. 55) TC CALCULATE THE BESSEL FUNCTION RATIO. THE P 
C AM: 0 COEFFICIENTS USED IN THi. EXPANSIONS ARt CALCULATED USING THE 
C SFfilES SUBROUTINE. 

C 

GO TO 20 

4 C X = R I- A L ( C A R G ) 

Y = -y 
DO 2 J = 1 , 2 
M = M + J 
XI = f LHAT ( k ' ) 

FM = C>PLX(XI,O.C) 

CALL ScRIES 
X P SUM ( J ) = P SUM 

2 XUSUM(J) = OSUM 

DO 3 K = 1 , 2 

A = X + (0.75 - FLOAT (K) ) * P I 

XX = XPSUM(K) * CCS ( A ) - XOSUM(K) » SIN(A) 

YY = (XPSUM(K) * SIMA) + XOSU’MK) » COS(A)) * Z 

3 X J ( K ) = XX + YY 

RATIO = X J ( 2 ) / X J ( 1 ) 

JR AT 10 = RATIO 
M = - 1 
GO TC JO 

2 C CAI I ZEEbES ( JO, O.CARC) 

CALL ZcFhLS(J2, 2, CARD) 

JR AT If. = J2/J0 
3 C CONTINUE 

RETURN 
END 
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tIBFIC SUBPUI I.IST 

COMPLEX FUNCTION PFI(J) 

THIS PROGRAM CALCULATES A FUNCTION PHI WHICH IS AN 
ATTENUATION PARAMETER. 

C 0 M M 0 \ / b E 0 K T Y / D ( 1C) ,XL( 10) , V (1 0 ) ,DVI)P(10) 

C OMM DM /ThER MO/GAMMA ,P ( 10 ) , RHO ( 1 0 ) , V I S ( 1 0 ) , A«PRES 
COMMOM/XFPEU/FR EC 
COMPLEX POLY , PH I , JR AT 10 

SN =6. 2M3185A«FRE(J*C( J ) **2*RH0( J) / VI S ( J ) 

PHI =6. si HyA*F RFC* SORT (RHO ( J )/ ( GAMMA * AMPRE S ) ) * 

P CS' RT ( 1. /JRATIOI CSCKTI CMPLX(0.0,-( SNI/4.0) ) ) ) )* 

P CSQRT ( GAMMA/POLY ( J ) ) 

RETURN 
END 


IPTIC SUBPLY LIST 

COMPLFa FUNCTION POIY(J) 

THIS ROUTINE CALCULATES A POLYTPOPIC COEFFICIENT 
( LABEL ED 'N' IN PERCH AND TIDJEMAN) 

COMMON/ GtCMTY /D ( 1C) ,XL( 1 0 ) , V ( 1 0 ) , D VI) P ( 1 0 ) 

COMMON / ThER w O/G AM "A ,P ( 1 C ) , RHIJ ( 10 ) . V I S ( 1 0 ) , AMPRF S 

COMMOU/XFRr'i.'/FREC 

COMPLEX JRAT IO.CARG.POLY 

SN = A.2fc31354*FREQ*D(J)**2*RIIO(J)/VIS(J) 

ARG = P(J)«SN/4.r 
C A R G = C S U R T ( CMP LX (U .0 ,-ARG ) ) 

PniY = 1.0/(l.G+( ( CAMMA- 1 .0 ) /GAMMA ) * JRATI OICARG ) ) 

R F TURN 
END 



$ IBFTC FCSINH LIST 

COMPLEX FUNCTION CSINH(ARG) 
COmdlFX ARGtCSINh 

C S I N H = (CFXP(ARG)-Cr. XP(-ARG) )/2.0 

RETURN 

END 


$IBFrC FCCOSF LIST 

COMPLEX FUNCTION CCCSH(ARG) 

COMPLEX ARC » CCO SF 

CCnSF = (CLXPIARC )+CFXP(-ARO) )/2.0 

RETURN 

END 
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